Stabilizing Grafted Colloids in a Polymer Melt: 
Favorable Enthalpic Interactions^ 



Itamar Borukhov* and Ludwik Leibler 

Unite Mixte de Recherche CNRS/Elf Atochem (UMR 167), 
95 rue Danton, B.P. 108, 92303 Levallois-Perret Cedex, France 
(February 6, 2008) 

The interactions between spherical colloids covered with end-grafted polymers (brushes) im- 
mersed in a polymer melt are studied theoretically. We show that attractive enthalpic interactions 
between the two polymer species, characterized by a negative Flory parameter (x < 0), can stabilize 
the colloidal dispersion. The stabilizing mechanism is a result of a structural change from a dry 
brush, where the melt chains do not penetrate deeply into the brush, to a wet brush, where the free 
chains penetrate the brush and force the grafted chains to extend into the melt. 
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Introduction: Mixtures of colloids and polymers pro- 
vide an important class of advanced materials with a wide 
range of mechanical, adhesive and optical properties. In 
many applications polymer chains are end-grafted to the 
colloid surface forming what is usually referred to as a 
polymer brush When the grafted particles are im- 

mersed in solution steric interactions acting between the 
grafted chains stabilize the colloidal suspension against 
attractive van der Waals interactions. Unfortunately, 
when the embedding medium is a polymer melt this 
mechanism does not always work. Due to subtle entropic 
effects the melt chains are expelled from the grafted layer 
even if they are chemically identical to the grafted chains 
and the inter-colloidal interactions become attrac- 
tive |9| . This attraction limits our ability to produce new 
materials where the contrast between the colloidal and 
polymeric properties might have advantageous implica- 
tions. The aim of this letter is to put forward a simple 
mechanism of stabilizing the grafted colloids in the sur- 
rounding polymer melt by means of a favorable enthalpic 
attraction between the grafted and free chains. A number 
of such systems exists, e.g., poly (methyl methacrylate) 
(PMMA) and polyvinyl chloride) (PVC), poly(styrene) 
(PS) and poly(2,6-dimethyl-l,4-phenylene oxide) (PXE), 
etc. In practice, these specific interactions have been 
used, for example, in the design of impact modifiers. 
Here, we quantify this effect and provide simple criteria 
for predicting the inter-colloidal interactions. 

Colloidal systems are not the only example where en- 
thalpic interactions compete with entropy in determining 
the structure of the material. Similar behavior can be 
also observed in mixtures of homopolymers and diblock 
copolymers JlOfl . Block copolymer chains are formed from 
two (or more) different polymers that are tethered to- 



gether at a single point. For example, an AB diblock 
copolymer consists of a sequence of A monomers followed 
by a sequence of B monomers. Depending on the incom- 
patibility of the two blocks these systems have a tendency 
to segregate into A- rich and B-rich mesophases O] . The 
resulting phases include lamellar structures of alternat- 
ing A-rich and B-rich stripes, hexagonal phases of cylin- 
drical domains, cubic phases of spherical domains, etc. 
Triblock copolymers (e.g. ABA or ABC) exhibit an even 
larger variety of phases |lT|] . When homopolymers are 
mixed with copolymers the mixture may exhibit very in- 
teresting physical properties (e.g. mechanical, optical, 
etc.) and can be used as, e.g. super tough materials or 
photonic devices |]l2| . Many of these properties depend 
crucially on the special meso-structure of the material 
and we would like to be able to determine under what 
conditions homopolymers can be added without destabi- 
lizing the structure. 

Our model system consists of two grafted spheres im- 
mersed in a polymer melt (Fig. |l]a). Typically, the radii 
R of the colloids is of the order of microns while the at- 
tractive interactions appear at separations w of order of 
1 to 10 nanometers. Such a separation of length scales 
allows us to study first a simpler case of two flat surfaces 
at a distance w from each other (Fig. [I}d) and only then 
include curvature effects. Each surface carries a dense 
layer of a /a 2 end-grafted chains per unit area, a being 
the effective size of a single monomer. The surfaces are 
immersed in a melt of mobile polymer chains, the indices 
of polymerization of the grafted and free chains denoted 
N and P, respectively. In the present study we restrict 
the discussion to the limit of long free chains P > N and 
to high grafting densities, while the more general case will 
be discussed elsewhere O]. For simplicity we assume 
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that the different monomer species are all of the same 
size. We further assume that the different monomers 
have the same volume a 3 , so that the incompressibility 
of the polymer melt can be expressed as 4>n + 4>p = 1 
where 4>n(x) and 4>p(x) are the local volume fractions of 
grafted and free chains, respectively. The enthalpic in- 
teraction between the different monomer species (in units 
of the thermal energy ksT) is characterized by the Flory 
interaction parameter \. 

Parabolic Approximation: In order to demonstrate 
the effect of a negative \ parameter we implement a 
variation of the parabolic approximation || , valid in the 
strong stretching limit of the grafted chains. In this limit, 
the mean- field potential felt by a single chain \x^{x) is 
a quadratic function of the distance x from the nearest 
surface fhi]] . Taking into account the symmetry around 
the mid-plane one can write 

, x_J A-Bx 2 Q<x<w/2 m 
PUNK*) - | A-B{w-x) 2 w/2<x<w (1 > 

where B — 3ir 2 /8N 2 a 2 , A(w) is an integration constant 
to be determined later and f3 — 1/fcsT. This potential 
has to be equal to the sum of the different interactions 
of the grafted chain with its environment. The latter can 
be obtained from the excess free energy associated with 
the free chains Vp = (1/a 3 ) J dxvp(x) where 
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The first term is the elastic energy of the free polymer 
chains, the second term is their translational entropy and 
the last term is the enthalpic interaction between the 
grafted and free chains. 

Taking into account the incompressibility constraint 
the mean field potential can be obtained from the varia- 
tion of the interaction energy (1n(x) — —SVp/S4>N which 
yields 
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Please note that due to the incompressibility of the poly- 
mer melt the last term has the same form as an effective 
excluded volume contribution with u e g ~ 2|x|a 3 . The 
above equation combined with eq. |l| defines a differen- 
tial equation for 4>n{x) in the interval < x < w, to 
be solved numerically. Since there is one free parameter, 
A(w), three constraints are required. These are provided 
by mass conservation, Naa — J dx(f>N(x), and by two 
boundary conditions <p' N = at x — and <p' N = at 
x = w/2. The former reflects our choice of a neutral sur- 
face that is equally indifferent to the two polymer species, 



while the latter reflects the symmetry of the problem. 
The free energy of interaction per unit area can then be 
calculated from 



F P ai{w) = 
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The above expression results from a transformation from 
single chain trajectories to an average volume fraction 
4>n{x). The first term is the standard free energy of an 
order parameter in a local field hn{x) while the last con- 
tribution is a reminiscence of the single chain partition 
function j|. 

The effect of a negative \ parameter is demonstrated 
in Fig. ||a where the excess free energy of interaction 
AFpol(w) = F po i(w) — Fp i(oo) is plotted as function 
of the inter-surface separation w. When the two poly- 
mer species are chemically the same (namely, \ — 0) 
an attractive minimum appears at short distances. This 
attractive minimum disappears in the presence of at- 
tractive enthalpic interactions between the two species 
(x = —0.15 in the figure (l6j). This behavior is quite 
general and below we provide criteria for the transition 
between attractive and repulsive behavior. 

The interaction energy of two colloids can now be es- 
timated using the Derjaguin approximation []15[ 



AU R (w) = ttR / dw'AFp ol (w') 



(5) 



Typical inter-particle interaction energies are plotted 
in Fig. ^b for colloids of radii R = 400a. Clearly, the 
inter-particle attraction disappears when enthalpic in- 
teractions are included. 

Discussion: Since the nature of the inter-surface in- 
teractions depends strongly on the shape of the two 
brushes when they begin to overlap, it is instructive to 
consider a single brush in contact with a polymer melt 
In Fig. ||a two volume fraction profiles of the 
grafted chains are plotted corresponding to \ = and 
X = —0.15. Since at high grafting densities the free poly- 
mers are too large to penetrate deep into the grafted 
layer, the grafted layer forms a so-called dry brush whose 
width h ~ Naa is determined by packing constraints. 

The free polymers only penetrate a narrow layer of 
width A, which is determined by the balance between the 
energy required to squeeze some free chains into the dense 
brush and the elastic response of the polymer chains to 
local variations in the density. These can be estimated by 
assuming a hyperbolic tangent profile for the volume frac- 
tion of the grafted chains 4>n{x) = (1— tanh[(a;— h)/X\)/2. 
The elastic contribution to the free energy per unit area 
is then ifTol 
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To first approximation the stretching energy can be es- 
timated by calculating the energy difference between a 
sharp box- like brush density and the smooth one |6|JlCt|: 



32 Na 4 



Minimizing the total free energy F to t = F c \ + F 1 
respect to A yields 
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(7) 
with 
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For the parameters of Fig. ||a one obtains A ~ 5.15a. The 
corresponding hyperbolic profile is indicated in the figure 
by a dotted line. The good agreement with the numeri- 
cal profile demonstrates that indeed F e i and F pen are the 
dominant terms in the free energy. 

The effect of enthalpic interactions with a negative x 
parameter on the structure of the brush is somewhat 
analogous to replacing a theta solvent with a good sol- 
vent where mixing of solvent molecules with the grafted 
polymer chains is encouraged [Q. The enthalpic contri- 
bution to the free energy of a dry brush can be easily 
estimated by 



X f X^ 

(3F X ~ / &X(j) N (x)(j)p(x) ~ 



(9) 



When \x\ is large enough F x exceeds F e \ and the penetra- 
tion length A results now from the competition between 
F x and F pon . Thus, for \ X \ > (tt 8 /48)^ 3 (a/N) 2 ^ the 
brush enters the enthalpic dry brush regime (see diagram 
in Fig. ||b) where the size of the brush is still h ~ Naa 
but the penetration length follows a different scaling law 
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As the enthalpic interactions become stronger, the pen- 
etration length A increases until A ~ h. This happens for 



(11) 



At this point the free chains penetrate all the way 
through the brush while the grafted chains further stretch 
into the polymer melt. The brush is now in the enthalpic 
wet brush regime (see Fig. |b). 

In this regime the enthalpic term dominates the 
potential pln{x) (eq. |^) and the density profile is 
parabolic, 4>n{x) = B(h 2 — x 2 )/2\x\ while h — 
(8/7r^) 1 / 3 (|x|cr) 1 / 3 A r a. Indeed, for the parameters of 
Fig. ||a with X = —0.15 one obtains h ~ 91a. 

When |x| is further increased beyond \x\ = 7r 2 /8<r, 
the grafted chains reach their maximal possible length 
/i max = Na and we enter the fully stretched brush regime 
(Fig. ^b). The parabolic approximation can be modified 
in order to take into account the finite extensibility of 
the chains |l7]] but this is beyond the scope of the cur- 
rent work. 



Comparison with van der Waals interactions: The re- 
pulsion between the grafted layers should be now com- 
pared with the van der Waals interactions. These in- 
teractions result from sharp variations in the dielectric 
properties across an interface separating two materials. 
When two surfaces of material (1) are separated at a 
distance I by media (2) the non-retarded van der Waals 
interaction energy can be written as []i~5"f 
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The Hamaker constant A121 depends on the static di- 
electric constants of the two media and their refractive 
indices. When the surfaces as well as the intervening me- 
dia are polymers the contrast is usually not very strong 
and the Hamaker constants are typically of order 10 -21 J 
(recall that lk B T ~ 4 x 10~ 21 J). 

When two dry brushes interact with each other the 
contrast at the brush/melt interfaces should be consid- 
ered in addition to the brush/surface interfaces. The 
former gives a stronger contribution since the separation 
between the two brushes I = w — 2Na is smaller than 
the bare inter-surface separation. In contrast, when two 
wet brushes interact with each other the brush/melt in- 
terface is no longer well defined because of the parabolic 
concentration profile which decreases rather moderately 
to zero (see Fig. ||a). Furthermore, since the free melt 
chains are present inside the brush the contrast is much 
smaller. This leaves only the van der Waals attraction 
between the bare surfaces which is too weak at the dis- 
tance imposed by the swelling of the brush. 

Conclusions: In this paper we have suggested a sim- 
ple mechanism for stabilizing grafted colloids in a poly- 
mer melt. We have shown that a negative x param- 
eter for the interaction between the grafted polymer 
chains and the mobile melt chains will induce swelling 
of the brush and consequently strong inter-brush repul- 
sion. This mechanism can have various implications for 
creating novel polymer-based materials. 
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FIG. 2. (a) Excess free energy of interaction per unit area AF po i between two fiat surfaces (in units of fcsT/a 2 ) and (b) 
interaction energy of two grafted colloids AUr (in units of ksT) as function of the reduced inter-surface distance w/a. The 
two solid curves correspond to \ = with grafting densities a — 0.02 and a = 0.05. The two dashed curves correspond to the 
same grafting densities with a negative enthalpic interaction \ = —0.15. The polymerization indices are N = P = 500. The 
radii of the colloids in (b) is R=400a. For clarity, the distances of closest approach (w — 2Naa) are indicated by vertical lines 
at w = 20a and w = 50a. 




FIG. 3. (a) Volume fraction profiles 4>n for x = (solid curve) and \ = —0.15 (dashed curve) as function of the distance 
from the surface x/a. The physical parameters are a = 0.05 and N = P = 500. The dotted line corresponds to the hyperbolic 
profile with A = 5.15. (b) Sequence of brush states as function of the negative Flory parameter \- See text for more details 
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